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Theoretical treatments of the BCS-BEC crossover need to provide as accurate as possible descrip- 
tions of the two regimes where the diluteness condition applies, either in terms of the constituent 
fermions (BCS limit) or of the composite bosons which form as bound-fermion pairs (BEC limit). 
This has to occur via a single fermionic theory that bridges across these two limiting representa- 
tions. In this paper, we set up successive improvements of the fermionic theory, that result into 
composite bosons described at the level of either the Bogoliubov or the Popov approximations for 
point-like bosons. This work bears on the recent experimental advances on the BCS-BEC crossover 
with trapped Fermi atoms, which show the need for accurate theoretical descriptions of BEC side 
of the crossover. 



PACS numbers: 03.75.Ss, 03.75.Hh, 05.30.Jp 



I. INTRODUCTION 

Recent experimental advancesi*^ with ultracold 
trapped Fermi atoms have produced systems of compos- 
ite bosons (dimers), where pairs of fermions bind due 
to their strong mutual attraction induced by a Fano- 
Feshbach resonance'^. Contrary to the case of point-hke 
bosons, for which the internal fermionic structure is im- 
material, for the composite bosons the internal fermionic 
structure is relevant to the extent that their binding 
energy is comparable with the energy and temperature 
scales involved in the experiments. Theoretical descrip- 
tions of the composite bosons should thus take into ac- 
count not only their overall bosonic structure (associated 
with their center-of-mass motion), but also their com- 
posite nature in terms of the degrees of freedom of the 
constituent fermions. 

Any sensible theoretical approach to the BCS-BEC 
crossover should rely on a single theory, that recovers 
controlled approximations on both sides of the crossover 
and provides a continuous evolution between them. In 
a previous paper*, the approach originally developed 
by Popov for a weakly-interacting (dilute) superfluid 
Fermi gas^ was improved by including the effects of the 
Bogoliubov- Anderson modc^ in the diagonal fermionic 
self-energy. With this generalization, it was also showni 
that, in the strong-coupling limit of the fermionic attrac- 
tion, the same theory is able to describe a dilute system of 
composite bosons within the Bogoliubov approximation^. 

Theoretical work on the condensates of dilute trapped 
Bose gases^iS has shown that an accurate description of 
the experimental data can be obtained by using the so- 
called Popov bosonic approximation^. This approxima- 
tion leads to improvements over the Bogoliubov approx- 
imation, in that it treats the densities of noncondensed 
and condensed atoms on the same footing, at least in 
the temperature range (below the critical temperature) 
where the two densities are comparable to each other. It 
has been shown that the use of the Popov bosonic ap- 
proximation is especially relevant to the thermodynam- 
ics of trapped Bose atomsS. It leads, in particular, to 



a downward shift of the critical temperature-^", which 
has been recently measured with great accuracyi-^. The 
Popov bosonic approximation was originally conceived 
for a homogeneous gas^; its extension to an inhomoge- 
ncous gas was described in detail in Ref. IT^ 

Purpose of the present paper is to devise a fermionic 
approximation for the BCS-BEC crossover in the broken- 
symmetry phase, which in the strong-coupling limit of 
the fermionic attraction recovers the Popov approxima- 
tion for the composite bosons. We shall arrive to this 
Popov approximation by successive improvements of the 
fermionic theory, starting initially from the BCS mean- 
field approximation and then including fluctuations over 
and above it at different levels of sophistication. The 
theory described already in Ref. (see also Ref. 1^ will 
be recovered as an intermediate step of this process. It 
will turn out that the Popov description of the compos- 
ite bosons on the strong-coupling side of the crossover 
corresponds to including pseudo-gap effects on top of the 
BCS approximation in the fermionic two-particle Green's 
function on the weak-coupling side of the crossover. Im- 
provements of the theory on the strong-coupling side, 
as to get a more accurate description of the scattering 
processes between the composite bosons, will be further 
considered along the lines of Ref. The issue to refine 
the theory on the weak-coupling side, as to include all 
possible contributions to an imperfect Fermi gaai^, will 
instead be the subject of future study. 

In the following, the theory of the BCS-BEC crossover 
will be considered for a homogeneous system. Extension 
to trapped gases may be done, in practice, by consider- 
ing a local-density approximation, whereby the fermionic 
chemical potential fj, is replaced whenever it occurs by the 
local quantity /x(r) = ^ — V{r) which includes the poten- 
tial V{r) at position r in the trap. This procedure has 
been already implemented at the BCS (mean-field) level 
in Ref. as well as with the inclusion of fluctuations 
in Ref. [13 (leading to the Bogoliubov description of the 
composite bosons). 

Our approach is built on a many-body Hamiltonian 
for fermions only, with an effective mutual attraction 
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represented by a point-contact interaction. It thus con- 
trasts an alternative approach based on a boson-fermion 
mode&i^^. It has, in fact, been shownSlthat (at least 
for the broad resonances currently used in most exper- 
iments) a point-contact interaction is by far sufficient 
to get an accurate description of the relevant scattering 
properties. 

The present paper deals with several kinds of approx- 
imations for the BCS-BEC crossover at a formal level. 
In this respect, alternative strategies will be proposed 
which emphasize the different constraints one has to deal 
with to obtain a satisfactory description of the composite 
bosons in terms of the constituent fermions. Detailed nu- 
merical calculations at the mean-field level and also with 
the inclusion of pairing fluctuations, which recover the 
Bogoliubov description of the composite bosons, have al- 
ready been presented in Ref. 0. Quantitative comparison 
among the different kinds of additional approximations 
discussed in the present paper awaits further numerical 
calculations, to be considered separately. 

The plan of the paper is as follows. In Section II, the 
steps leading to the Bogoliubov description of the com- 
posite bosons (discussed already in Refs. and IT^h are 
reproposed, aiming at emphasizing the hierarchy of ap- 
proximations for the fermions and the composite bosons. 
The role played in this context by the gap equation will 
be specifically discussed. In Section III, it will be shown 
how to modify the fermionic theory in order to describe 
the composite bosons within the Popov approximation. 
In both (Bogoliubov and Popov) cases, it will be also 
discussed how to enlarge the selection of the fermionic 
many-body diagrams, as to refine the description of the 
scattering processes between the composite bosons. Sec- 
tion IV gives our conclusions. For the sake of complete- 
ness, the Bogoliubov and Popov descriptions for point- 
like bosons will be briefly recalled in the Appendix. 



II. BOGOLIUBOV DESCRIPTION OF 
COMPOSITE BOSONS 

As anticipated in the Introduction, we will obtain the 
Popov approximation for the composite bosons by set- 
ting up a scheme of successive approximations for the 
constituent fermions. To this end, it will be useful to re- 
trace the essential steps which have previously led to the 
description of the composite bosons within the Bogoli- 
ubov approximation (cf. Refs.0 and lT3|) . By this proce- 
dure, the interplay between the approximations adopted 
for the constituent fermions (which emerge especially on 
the BCS side of the crossover) and the ensuing approxi- 
mations for the composite bosons (which are mostly rel- 
evant to the BEC side of the crossover) will become evi- 
dent. The material presented in this Section also serves 
for setting up the notation that will be needed in the fol- 
lowing Section, where the Popov approximation for the 
composite bosons will be explicitly introduced. 

As also mentioned in the Introduction, the present 



approach is based on a many-body Hamiltonian for 
fermions only, with the (effective) mutual attraction rep- 
resented by a point-contact potential vq S{r) where vq is 
a negative constant. This choice entails a suitable regu- 
larization in terms of a cutoff ko in wave- vector space, by 
taking the constant vq of the form (in three dimensions) 

_ _ 27r^ _ ^3 

where m is the fermion mass and a p the scattering length 
of the associated (fermionic) two-body problemi4. In 
this way, the relevant information about the two-body 
problem is fed into the many-body problem. By letting 
ko oo (and thus vq 0) eventually, the classification 
of the (fermionic) many-body diagrams gets considerably 
simplified in the normal phase^** as well as in the broken- 
symmetry phase^^, since only specific diagrammatic sub- 
structures survive in the limit. 



A. Fermions within the BCS approximation 

The BCS mean field represents the simplest approx- 
imation for the BCS-BEC crossover in the broken- 
symmetry phase. As such, it is often taken as a ref- 
erence approximation for this crossover, at least at low 
enough temperatures^. With this approximation, the 
fermionic single-particle Green's functions (in Nambu's 
notation) are obtained by solving the Dyson's equation 
in the formi: 

(?fcs(p) ^ -(7BCS(_p) ^ g^^p) ^ 5o(p)(-A)ei\cS(p) 

eiT'(p) = Gfi'^Hp) = go{pK-A)gf,^'{p) ■ (2) 

Here, p = (p,^^) is a four-vector with wave vector p 
and fermionic Matsubara frequency cug = (2s -I- l)7r//3 (s 
being an integer and /3 the inverse temperature), Goip) = 
(iujs — p^ /{2m)+^)~^ is the free-fermion propagator, and 
we set h = kB = ^ throughout. The isotropic (s-wave) 
BCS gap function A is obtained by solving the equation: 

'^—/^iE^ff^rt. (3) 

The regularization of the contact potential is provided in 
this context via the relation 

TO 1 f dp m 

4^ " W ^ J J2^^ ^ ' 

which defines the scattering length ap and is equivalent 
to Eq. Q. One is led in this way to introduce the di- 
mensionless coupling parameter {kpap)~^ (where kp is 
the Fermi wave vector related to the density n via n — 
kp/{37T'^)). This parameter ranges formally from — oo in 
weak coupling to -l-oo in strong coupling, although, in 
practice, the crossover occurs for {kp\ap\)^^ < 1. 

Besides being relevant per se to the BCS-BEC 
crossover, in the present context the BCS gap equation 
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(O serves to determine the bosonic chemical potential /is 
which enters the Bogoliubov approximation for the com- 
posite bosons in the strong-coupling limit of the fermionic 
attraction. This is because the BCS gap equation Q 
plays an analogous role to the Gross-Pitaevskii equation 
for point-like bosonsSi. For a homogeneous system, it 
thus fixes the relation between fis and the condensate 
density ng. 

To show this, we perform the sum over the Matsubara 
frequency in Eq. ||2J) and eliminate Vq using Eq. to 
obtain: 



dp 

(2^ 



tanh(/5£;(p)/2) 



m 
7^ 



(5) 



where £:(p) = ^^(p)2 + A2 with ^(p) = p'^/{2m) - 
fi is the BCS quasiparticle dispersion. In the strong- 
coupling (A ^ limit^, one can approximate the 
hyperbolic tangent in the expression 101 by unity and 
expand E{p)^^ to first order in (A/|/i|)2. The resulting 
integrals can be done analytically, yielding^: 



4|mI 



fJ-B 



(6) 



Here, the formal relation /is — 2/x -I- eg between the 
bosonic and fermionic chemical potentials'^ has been 
used, where ep = (rna^)~^ is the binding energy of the 
fermionic two-body problem. 

The condensate density no for composite bosons is fur- 
ther identified as jaj^, where in the BCS approximation 



IPJ = \ — 




(7) 
(8) 



being the (normalized) wave function of the fermionic 
two-body problem. Using the Green's function one 
then obtains in the strong-coupling limilji^: 



a ~ A 



im?ap 



(9) 



Correspondingly, the gap equation (0 reduces to the 
form 



/ Anas \ 
V riiB } 



(10) 



where ms = 2to is the mass of the composite bosons 
and as — 'i.ap is the value of the scattering length for 
the residual boson-boson interaction within the present 
approximation. [We shall discuss below how improved 
descriptions of the boson-boson scattering result into 
smaller values of the ratio aB/ap-] A relation between 
/is and no similar to H1U|) holds for point-like bosons 
within the Bogoliubov approximation. Its enforcement 
on the Bogoliubov quasiparticle dispersion leads to a 



gapless spectrum, as required on general ground by the 
Hugenholtz-Pines theorem^. 

The strong modification of the fermionic chemical po- 
tential /X when passing from the weak-coupling BCS limit 
(where /i ~ tp, tp = k'^p / {'^ni) being the Fermi energy) 
to the strong-coupling BEC limit (where /i ~ — eo/2), 
quite generally, results by supplementing the BCS gap 
equation ((Sj with the density equation: 



n = 2 



dp 1 



(27r)3 [3 



(11) 



where 7^ = 0"'" and the factor of 2 accounts for the two 
equally populated spin components. In particular, by 
expanding the BCS single-particle Green's function ((SJ to 
the leading significant order in A, in the strong-coupling 
limit one obtains: 



dp 1 



(27r)3 13 



(12) 

The first term on the right-hand side of Eq. (|12|l vanishes 
in the strong-coupling {f3fi — > — oo) limit. In the second 
term we instead approximate 



dp 1 



(27r)3 (3 



Y.g^{pfg,{~p) 



m? ap 



so that Eq. H12|) reduces to 



n . 9 m ap , 
_^A^__^H2 



(13) 



(14) 



according to Eq. 0. No depletion of the condensate, 
therefore, occurs within the BCS mean-field approxima- 
tion. Physically, this corresponds to temperatures much 
smaller than the critical temperature and to a weak resid- 
ual boson-boson interaction. To describe the composite 
bosons at temperatures comparable with the critical tem- 
perature, it is then required to include fiuctuation correc- 
tions beyond BCS in the fermionic single-particle Green's 
functions, as it will be discussed in subsection IIC. 



B. Composite-boson propagators within the 
Bogoliubov scheme 

With the fermions described by the BCS approxi- 
mation, one can construct the ladder propagators for 
fermion pairs as described in ReflU One obtains: 



rii((?) Ti2{q) 
r2i(g) T22iq) 



Aiq)Ai-q) - B{qY 

■ A{-q) B{q) 
B{q) A{q) 



(15) 



Here, q = (q, fli^) is a four-vector with wave vector q and 
bosonic Matsubara frequency ^li, = 2ttv/I3 {v integer). 
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FIG. 1: (a) Graphical representation of the normal (left) and 
anomalous (right) rungs entering the ladder propagators of 
Eq. 1151 . (b) Normal and (c) anomalous self-energies for the 
composite bosons resulting from Egs. 12411 and 12511 . respec- 
tively. Thick lines stand for the single-particle BCS Green's 
functions @, light lines for the free-fermion propagator Qo, 
and black squares for (minus) the gap function A. 



and 



1 



dp 1 v;^^ 



BCS 



'Z)aiT'(-p) (16) 



contain the rungs depicted in Fig. 1(a). Only these two 
rungs need, in fact, be considered owing to the regulariza- 
tion ^ of the fermionic potential. In Ref. it was also 
shown that, in the strong-coupling limit of the fermionic 
attraction, the ladder propagators 1)15(1 reduce to: 

rii(g) T22{-q) 



Ti2{q) = T2i[q) ~ 
where 




(17) 
(18) 

(19) 



is the Bogoliubov quasiparticle dispersion with fiB given 
by Eq. (HHJ. The results 1(13 and (HHJ bear strong 
resemblance, respectively, with the normal {G'b) and 
anomalous {G^) noncondensate bosonic Green's func- 
tions within the Bogoliubov approximation^. The com- 
parison gives: 



rii(g) = -{87r/m^aF)G's{q) 
Ti2{q) - i8n/m^aF)Gl\q) 



(20) 



where the sign difference between the two expressions was 
accounted for in Ref. 13. 

The form l|15|l of the ladder propagators is obtained 
by solving the Bethe-Salpeter equation for the fermionic 
two-particle Green's function in the particle-particle 
channel. It thus holds quite generally for any value of 
the fermionic coupling {kpap)~^ . The rungs in the ex- 
pressions (I16II contain the BCS single-particle Green's 
functions which are the self- consistent solutions of the 
Dyson's equation ((2Jl. This is in accordance with the pre- 
scriptions for the approximation to be "conserving"i^Si2i. 
In addition, the condition ^(0) — -B(O) ~ 0, which is re- 
quired for the ladder propagators ((T^ to be "gapless" at 
g = 0, is equivalent to the BCS gap equation © for any 
coupling. The equivalence can be readily proved via the 
identity 

e,TS(p) = A [erT^(p)afP(-p) + erT^(p)aiT^(-p)] 

(21) 

which holds for the BCS Green's functions ||2Jl. It im- 
plies, in particular, that the approximation H17|) and l|18|) 
for the propagators of the composite bosons is conserving 
and gapless. This result is consistent with a general prop- 
erty proved in Ref. |2^ according to which a given con- 
serving approximation for the constituent fermions also 
results into a gapless approximation for the composite 
bosons (provided the Baym-Kadanoff prescriptions are 
satisfied also in the two-particle channel). 

In strong coupling the above condition of self- 
consistency can, however, be relaxed since the system 
of composite bosons becomes dilute. The presence of the 
small "gas parameter" n^J'^UB = rkFaF{6Tr'^)~^^^ (where 
nB = n/2 is the bosonic density and the ratio r = as/ap 
is of the order unity) suffices, in fact, to select the rele- 
vant diagrammatic approximations for a dilute Bose gas 
(barring the temperature regions close to the critical tem- 
perature and to zero temperature)^. 

In this case, the forms (|17|l and (|18|) for the propa- 
gators could be obtained by solving a Dyson's equation 
directly for the composite bosons, in an analogous way 
to Eq. (|A1I) for point- like bosons, with the ladder propa- 
gator To{q) for the normal phase playing the role of the 
free-boson propagator. At arbitrary coupling this ladder 
propagator takes the forrr^-: 



ro(g) 



r m 



dp 



tanh(/3C(p) /2) + tanh(/3C(p - q) /2) m 



(22) 

-1 



2(e(p) + e(p - q) - 1^.) P 

In particular, in the strong-coupling limit it reduces to 
87r 1 



ro('z) 



m?aF i^u — q^/(4TO) + /is 



(23) 



which coincides (apart again from the overall factor 
— 87r/(m^ai?)) with the free-boson propag ator G'-°\q) of 
the Appendix. 
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The corresponding normal and anomalous self-energies 
for the composite bosons were identified in Appendix B 
of Ref. IT^ by expanding the expressions of the rungs 
which result after the sum over the Matsubara fre- 
quencies is performed, to the leading significant order 
in A. For the following purposes, it is convenient to 
rephrase that analysis by expanding directly the BCS 
single-particle Green's functions Q entering the rungs 
(|16|l to the leading significant order in A. For these rungs 
one thus obtains: 

y-Q«{p + q)Q„(-p)-i^'I^^Y.9M^S«{-rf (24) 

and 

-A^/(^^E^o(p)^ao(-p)^ (25) 

with the neglect of the g— dependence whenever irrele- 
vant. Note the presence of a factor 2 on the right-hand 
side of Eq. (|24l) . This factor originates from the two 
distinct contributions to the normal self-energy for the 
composite bosons depicted in Fig. 1(b), which result upon 
dressing the upper and lower fermionic lines, respectively. 
The anomalous self-energy of Fig. 1(c), on the contrary, 
contains only one contribution (also with a sign difference 
with respect to the normal self-energy). 

In the expressions (|24|) and H25|l . the quantity 

^2(0, 0,0,0) ^ ^ g,{pfg,{-pf (26) 

can be interpreted (apart from an overall factor, see 
Eq. H28|l below) as an effective boson-boson interaction 
where all incoming and outgoing four-vectors vanish. 
This quantity thus represents a degenerate form of the ef- 
fective boson-boson interaction for arbitrary four-vectors 
(depicted in Fig. 2(a)) 



"2(^1, 92, 93, 94) = 




X Qo{-p)Qo{p + q2)Qo{-p + 91 - qi)Qo[p + 94) (27) 

with ch -l- 92 = 93 + 94, which was introduced in Refs. I25I 
andjlj- Note that the expression (|27|l contains the free- 
fermion propagator t/o of the normal phase (with an ap- 
propriate value of the fermionic chemical potential - see 
below). This is consistent with the fact that also for 
point-like bosons (cf. the Appendix) the effective boson- 
boson interaction (when associated with a t-matrix) is 
considered to be mildly dependent on temperature and 
is correspondingly evaluated in the normal phase^. In 



P+q4 




+ — 



FIG. 2: (a) Graphical representation of the effective boson- 
boson interaction 112 of Eq. 1271 . (b) First additional terms 
contributed by the t-matrix ts for the composite bosons given 
by Eq. 1311 . Like in Fig. 1, light lines stand for the free- 
fermion propagator, while broken lines stand for the fermionic 
interaction potential. For simplicity, spin labels are not shown 
explicitly. 

the strong-coupling limit, one obtains for Eq. (|26|l : 

Apart from the overall factor (m^ai?/(87r))^ (which is 
needed to compensate the factors ~8tt/ {ni^ap) originat- 
ing from the expression H23|) associated with the free- 
boson propagator), the result H28|) is consistent with the 
residual boson-boson interaction entering Eq. H10|l and 
yields again the value 2 for the ratio r = as/aF- 

The above considerations suggest us how to improve on 
the Bogoliubov approximation for the composite bosons, 
by including, in particular, the diagrammatic contribu- 
tions which have been shown in Ref. to decrease the 
ratio r from the value 2 to about 0.75 in the zero-density 
limit^^. The price one has to pay for this improvement 
is to give up the self-consistency of the fermionic single- 
particle Green's functions, which characterizes the ladder 
propagators ifTCjl . 

To this end, we first obtain the ladder propagators 
Tb (9) approximately for any value of the fermionic cou- 
pling by adopting a Dyson's equation of the type IjAip . 
where now: 

(i) The propagators Gsiq) of point- like bosons are re- 
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placed by the ladder propagators ^siq) of composite 
bosons; 

(ii) In the inverse ljA2p of the free-boson propagator 
G^^(q), the expression (^i^u — + fJ-B^ is replaced 
by ro(q')"^ as given by Eq. if^ : 

(iii) The bosonic self-energy IjASp within the Bogoliubov 
approximation is replaced by 



-2u2{0,q,0,q) U2{0,0,~q,q) 
U2{0,0,-q,q) -2u2{0,q,0,q) 



(29) 

Note that the g— dependence of the diagrams of Fig. 1 
has been retained in the expression (|29|l . In the strong- 
coupling limit, however, this q— dependence can be ig- 
nored, and the diagonal and off-diagonal components of 
the self-energy (|29ll reduce to the expressions given in 
Eas. H24ll and (|25|) . in the order. The ladder propagators 
Tsiq) obtained with the above prescriptions (i)-(iii) re- 
duce correspondingly to the strong-coupling expressions 
(fT7| and 

Similarly to point-like bosons^ , the (bare) boson-boson 
interaction entering Eq. (|29|l can be conveniently evalu- 
ated in the normal phase, even though the correspond- 
ing self-energy refers to the broken-symmetry phase. 
This is because in the strong-coupling limit the fermionic 
chemical potential /i represents the largest energy scale 
in the problem-i, so that the precise value of the tem- 
perature T is irrelevant insofar as T ^ 

To guarantee that the ladder propagators Tsiq) are 
gapless at g = for any value of the fermionic coupling, 
the condition 

ro(g = 0)-^ - Sii(g = 0) - Si2(g = o) = O (30) 

needs to be satisfied. In the present context, this con- 
dition plays the role of the Hugenholtz-Pines theorem 
for point-like bosons!.. It can be shown that the condi- 
tion (|30|l reduces to the BCS gap equation (|3J| both in 
strong coupling (where it is equivalent to the result 
and in weak coupling. Accordingly, one expects it to ap- 
proximate the gap equation Q reasonably well for all 
couplings. 

The Bogoliubov approximation for the composite 
bosons can be improved at this point by considering the 
t-matrix tsiqi, (72, 93, (74) for the composite bosons intro- 
duced in Ref. 14 and obtained by solving the integral 
equation: 

^^(gi, 92, 93, 94) = ^2(91,92,93,94) 
~y-^^^^X!"2(9i,92,95,9i +92 - 95)ro(95) 

xro(9i +92 " 95)^5(91 +92 -95,95,93,94) • (31) 

The first few terms obtained by iteration of Eq. H31(l 
are represented in Fig. 2(b) in terms of the constituent 
fermions. These diagrams have been shown in Ref. ^3 
to correct the ratio aB/ap from the value 2 (obtained 



with the Born approximation for the boson-boson scat- 
tering) to about 0.75. The analysis made in Ref. |3 
also suggests that inclusion of these diagrams should be- 
come immaterial when approaching the crossover region 
(kpaF)~^ ~ 0, since these diagrams correspond to high- 
order fermionic pairing-fiuctuation processes^!? . To cor- 
rect for the value of the boson-boson scattering length 
entering the ladder propagators H17|) and 118|l in strong 
coupling, it is thus sufficient to replace U2 (91, 92, 93, 94) 
in Eq. by the t-matrix ^^(gi, 52, 93, 94) of Eq. J^J). 
Correspondingly, the condition Ij^UI) must be imposed for 
the composite bosons to be gapless. In strong coupling, 
this condition reduces to the result (|10|l with the modified 
value of a_B- 

The integral equation (|31|l was originally introduced in 
Ref. for the normal phase (that is, above the critical 
temperature Tc). Extension of that equation below Tj, 
would introduce infrared divergences owing to the pres- 
ence of the Goldstone mode^^. For point-like bosons, 
PopoAi^ has avoided these problems by considering the 
t-matrix associated with the two-body problem (that is, 
in the zero-density limit). In the present context, we 
may proceed along similar lines while preserving the com- 
posite nature of the bosons. To this end, we may pre- 
liminarly solve the integral equation H31|l at finite den- 
sity in the normal phase. For a given value of the cou- 
pling [kpap)^^ , the temperature is correspondingly kept 
(slightly) above Tc, where Tc and the chemical poten- 
tial /i are calculated according to the fermionic t-matrix 
approximation of Ref. 32. This value of /i is eventually 
inserted in Eq. (|31|l . 

What is still missing for a complete treatment of the 
Bogoliubov approximation for the composite bosons is to 
account for the depiction of the condensate, which is ab- 
sent in the BCS result (|14|l . This depletion occurs when 
approaching the crossover region from strong coupling 
and/or when the temperature is increased toward Tc. 



C. Fermions within BCS plus pairing fluctuations: 
Inclusion of the depletion of the condensate 

To find a relation between the condensate (no) and 
noncondensate (n') densities of the composite bosons, 
yet preserving the relation (|10l) characteristic of the Bo- 
goliubov theory, it was shown in Ref. 4 that the BCS 
density equation Hll|l needs to be modified as to include 
pairing-fiuctuation corrections, while preserving the BCS 
form (O of the gap equation. To this end, the fermionic 
Dyson's equation Q is replaced by 

Qii{p) = -Q22{-p) = Qa{p) + Qo{p)[^ii{p)Qii{p) 

+ Si2(p)a2i(p)] (32) 

Ql2{p) = Q2l{p) ^ Qq{p)[^Ii{p)Qi2{p) +'^12{P)Q22{P)] , 



where 

Sll(p) = -S22(-P) 



(9 -P) 



7 



q-p 




£„(P) 



P 



FIG. 3: Graphical representation of the fermionic self-energies 
1331 1 that include pairing fluctuations. The shaded box corre- 
sponds to the normal ladder propagator of Eq. 1151 . Conven- 
tions are like in Figs.l and 2. 



5]i2(p) = T,2i{p) = - A 



(33) 



with Fii given by Eq. (|15|l for arbitrary coupling. The 
fermionic self-energies (|33f) correspond to the diagrams 
of Fig. 3. The function Qff^ entering Eq. |j3Sl) has the 
functional form of the BCS single-particle Green's func- 
tion (0) albeit with modified values of A and /i, different 



from those obtained (at any given coupling) within the 
BCS approximation of subsection IIA. 

As anticipated, the two (gap and density) equations 
which drive the BCS-BEC crossover are treated here on 
a different footing*. The gap equation is still taken of the 
form with the BCS anomalous Green's function Grf^ , 
while the density equation now contains the dressed nor- 
mal Green's function Gii of Eq. and reads: 



dp 

(2^3 p 



(34) 



in the place of Eq. ifTTjl . 

The required relation between no and n' for the com- 
posite bosons can be obtained by expanding Qn in 
Eq. H34|) up to second order in the self-energies via the 
Dyson's equation (|32ll : 



Qii{p) = Qq{p) + 5o(p)Sii(p)t/o(p) 
+ ^o(p)Sii(p)eo(rtSii(p)t/o(p) 
- Go{p)^i2{p)Qo{-p)^2i{p)Qo{p) 



(35) 



which corresponds to an expansion in inverse powers of 
\^\. With the expressions H33|) for the self-energies, the 
density equation H34|l then reduces to: 



n 
2 



(36) 



Note that the g— dependence of the single-particle 
Green's function entering the expression (|33|l for En 
has been neglected. Correspondingly, Fn is taken of 
the polar form (|17|l which is valid in the strong-coupling 
(/3^ — > —oo) limit. The first term on the right-hand side 
of Eq. H36|l vanishes in this limit. In the second term, we 
write 



(2^)3 p . 



Stt 



(2^ 



m^ap 



(37) 



owing to the first of Eas. (|2()|l and to the standard def- 
inition of the noncondensate density n' for point-like 
bosons^. Using further the result H13|l which is also valid 
in strong coupling, as well as the related result 



647r 



(38) 



the expression H36() reduces eventually to: 



n 
2 



37rn' a|i) 



A' 



2 m ap 
Stt 



where 



<C 1 in the dilute limit of interest. 



(39) 



There remains to show that the second term on the 
right-hand side of Eq. (|39|l represents the condensate den- 
sity riQ even when pairing fluctuations are included, in 
agreement with the result that was proved within the 
BCS approximation. To this end, in strong coupling we 
approximate the diagonal self-energy (j33|l in the form: 



Sii(p) 



1 ^e^a..rn(,)) GU-v) 



- Is? 



;27r)3 /3 
1 



where we have used the result (|37|l for n' and set 



(40) 



8 



^0 



ap 



(41) 



With this expression for the self-energy, the fermionic 
single-particle Green's functions solutions of the Dyson's 
equation H32() are given by: 

Gnip) - K + e(p)-S22b)]/{[«t^s-e(p)-Sii(p)] 
iuJs + C(p) 



-7(P)^-A2 



(42) 



and 



GMp) = Si2(p)/{K-e(p)-Sii(p)] 

X K +?(p) - I]22(P)] - Si2(p)I]2l(p)} 

A 



(43) 



+ 7(P)^ 
with the notation 

7(p)2^e(p)' + A' + 2A2. (44) 
Using at this point the expression H43I) in the definition 
dp 1 



(45) 



which generalizes Eq. iQ to the present context, we ob- 
tain: 



a 



A f dp <f>{p) 
2 J (27r)3 7(p) 



A 



m?ap 



(46) 



This expression coincides with the BCS result lO to the 
leading order in A/|yLt|. 

Entering eventually the result H4t)|) into Eq. (|39|l and 
identifying again |ap with the condensate density ng for 
the composite bosons, we get the desired relationship 



n 
2 



no 



(47) 



which enables us to determine the coupling and temper- 
ature dependence of rip in terms of that of n' within the 
Bogoliubov approximation. Recall that the Bogoliubov 
result (|10() for the bosonic chemical still holds, since we 
have kept the BCS gap equation ((2Jl unchanged. 

As a final remark, we mention that the above results 
can also be obtained when considering the refinements 
discussed in subsection IIB to get a value smaller than 2 
for the ratio as/flF- 



III. POPOV DESCRIPTION OF COMPOSITE 
BOSONS 

As emphasized in the previous Section, the Bogoli- 
ubov approximation rests on the relation (|10|l between 



Hb and no. For point-like bosons this relation results 
from the Gross-Pitaevskii equation (in the homogeneous 
case) , while for composite bosons it is obtained from the 
BCS gap equation as shown in subsection IIA. For 
this reason, in subsection IIC the BCS anomalous Green's 
function Q^^^ was still utilized in the gap equation, even 
though the dressed normal Green's function Qn entered 
the density equation. 

In the Popov approximation for point-like bosons, on 
the other hand, the chemical potential is given by the 
relation ^lb — g{no + 2n') (cf. the Appendix), a result 
which can also be derived from a generalized version of 
the Gross-Pitaevskii equation^^. To get a similar relation 
for the composite bosons, it is thus clear that the BCS 
gap equation (|SJ) needs to be modified. 

To this end, we enter the expression H43|l for the 
anomalous fermionic single-particle Green's function 
(that contains pairing-fluctuation corrections beyond 
mean field) into the gap equation 



-Vo 



dp 

(2^3^ 



■Xp) 



(48) 



which generalizes Eq. (PJ- Considering specifically the 
strong-coupling limit and taking into account the regu- 
larization I^J for the contact potential, Eq. H48|l becomes: 



An a I 



dp 

(2^ 



tanh(^7(p)/2) m 
2^Kp^ ^ 



(49) 



which resembles Eq. ^ with E{p) replaced by 7(p) of 
Eq. (|44|l . In strong coupling, the hyperbolic tangent can 
again be approximated by unity, so that an expression 
similar to JBJ still results but now with A'^ replaced by 
A2 



2Ag: 



2A§) 



(50) 



Since the result (|46|l holds to the leadijng order in A/|/x| 
even when Q12 is dressed by pairing fluctuations, Eq. H5Ufl 
can be cast in the form 



\ ruB 



(no + 2n') 



(51) 



provided the relation (|37|l for n' still holds even within 
the present Popov approximation (as it will be verified 
below). Note that Eq. H51|l implies again the Born- 
approximation result as = 2a p- 

The next step to obtain the Popov approximation for 
the composite bosons is to identify additional diagrams in 
the fermionic particle-particle channel, which reproduce 
in strong coupling the additional (diagonal) bosonic self- 
energy (|A6I) of the Appendix needed to obtain the Popov 
approximation for point-like bosons. Let us similarly call 
S^°'^((7)ii and YF^'^{q)22 the corresponding self-energies 
for the composite bosons. The diagram associated with 
the dressing of the upper fermionic line in the ladder 
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q-q-p 



p+q 




p+q 



-P 

Pop 
£.(q) 



11 



FIG. 4: Graphical representation of the bosonic self-energy 
Eg°''(g)ii of Eq. ib'21 . obtained by dressing the upper 
fermionic line in the particle-particle channel. An analogous 
dressing can be done for the lower fermionic line. Conventions 
are like in Fig. 3. 



propagator is depicted in Fig. 4 (a similar diagram can 
be drawn for the dressing of the lower fermionic line) . It 



corresponds to the expression: 

y<Giiip + qfGii{-p)gii{q'-q-p)Tn{q') • (52) 

[The factor of 2 in Eq. (|52() accounts for the dressing of 
the the lower fermionic line, although this is not shown 
explicitly in Fig. 4.] In particular, in strong coupling we 
may approximate the diagram of Fig. 4 as follows: 



nPop/ 



\q)ii = ^B°^i-q) 



(27r)3 (3 



(53) 



where Fn is taken of the polar form H17|) as in Eq. H36|l. 

The Popov propagators for the composite bosons are 
now obtained in terms of the corresponding Bogoliubov 
propagators H15|l and of the above bosonic self-energy 



E^°^(g)ii , similarly to what is done in Eq. (|A7(I for point- 



like bosons. The result is: 



rP°P(g) rP°p(<z) 
rrr(g) r^r('z) 



/A(-g)-I]P°P(-g)n 



[A{q) 



nPop 



{q)ii][A{-q) 



^Pop 



-'B 
B{q) 



B{q) 



A{q) 



r 



"(9)11 ) 

m 



where A{q) and B{q) are given by Eq. l(TH|. 
The propagators are gapless provided 



A{q = 0) 



^Pop 



(g = 0)ii-B(q = 0) = 



(55) 



This equation generalizes to the present context the con- 
dition A{q = 0) — B{q = 0) = for gapless Bogoliubov 
propagators. Note, however, that while the condition 
A{q = 0) — B{q = 0) = is fully equivalent to the 
BCS gap equation for any coupling, its generaliza- 
tion (|55|l reduces to the gap equation (|48|l only in the 
strong- and weak-coupling limits. In strong coupling, 
this can be verified by writing 7(p)'^ = £'(p)^ + 2Aq in 
Eq. H49|) according to the definition H44|l . by expanding 
the resulting expression to first order in the small quan- 
tity (Ao/i?(p))^, and by recalling the expression (|^ . 
In weak coupling, on the other hand, the fermionic self- 
energy Eii(p) of Eq. (|33|l can be split into a constant 
part Sq and a remainder E{\(p), which depends explic- 
itly on p and whose magnitude is much smaller than A. 
By expanding the expression (|32|) for Q12 to the lowest 



significant order beyond BCS, one obtains in this way: 

^s+e(p) 



gi2(p)=^eiT'(p)-AI]fi(p) 



-ASfi(-p) 



iujs - g(p) 

[iUsf - i?(p)2 



(56) 



where the chemical potential has been shifted as /i ^ 
/i — So whenever it appears. Correspondingly, the gap 
equation H48|l becomes: 



1 

Wo 



dp tanh(/3£;(p)/2) 
(27r)3 2E{-p) 
dp 1 



+ 2 j ^i^S]n(p)eo(p)'eo(-p) = . (57) 

Note that the full Sii(p) has be restored in Eq. H57|l . since 
the last term on its right-hand side would vanish if En (p) 
were taken to be a constant, according to the approxima- 
tions vahd in the BCS hmit. With the aid of Eq. in 
the last term on the right-hand side of Eq. (|57|l one rec- 
ognizes the weak-coupling limit of the expression (|52|l for 
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the Popov bosonic self-energy. The weak-couphng hmit 
of the gapless condition H55|) is thus recovered from the 
gap equation (|48(l . Recall in this context that the need for 
including the constant shift Sq in the fermionic chemical 
potential entering the weak-coupling expressions of the 
rungs (|16|l was emphasized in Ref. 4 for all temperatures 
in the broken-symmetry phase. 

At arbitrary coupling, the broken-symmetry require- 
ment for the particle-particle propagator to be gapless 
can thus be satisfied by imposing the condition H55|l in 
the place of the gap equation H49|) . As shown above, this 
replacement does not affect quantitatively the weak- and 
strong-coupling limits of the theory. 

Caution must be exerted when implementing the form 
()52|l of the bosonic self-energy, where the bosonic propa- 
gator Fii appears. This propagator was originally meant 
to be taken of the form (|15|l . which is valid within the 
Bogoliubov approximation and to which the gapless con- 
dition A{q = 0) — B{q = 0) = applies. However, to use 
the alternative gapless condition H55|) for the Popov prop- 
agators (|54|l as well as for the bosonic self-energy (|52ll on 
which these propagators are built, Fn in Eq. (|52|) must 
be interpreted as given by Eq. H54|) itself. This procedure 
introduces a self-consistency in the bosonic calculation. 
By doing so, the above results are valid to the extent that 
Eq. I|37|) holds also within the Popov approximation. 

An alternative approach, which does not require reach- 
ing self-consistency in the bosonic expressions, would be 
to maintain the Bogoliubov form (|15|l for the propaga- 
tor Fii entering the bosonic self-energy (O with the 
associated gapless condition A{q = 0) — B{q = 0) = 0, 
while imposing the modified gapless condition (|55|l on 
the Popov propagators H54|l. 

The Popov results for point-like bosons (cf. the Ap- 
pendix) are fully recovered in strong coupling, since in 
this limit 



B{q) 



Stt 
m?ap 

Stt 
m?ap 

Stt 



2m f 



Snap 
no 



m 

Strap 
r 

TO 



-no 



(58) 



To get these expressions, Eqs.lQBJ and ((5^ have been uti- 
lized together with the generic definition /i^ = 2/x 4- eo of 
the bosonic chemical potential in terms of the fermionic 
one. Use of the specific form (|51f) for fis, as required by 
the condition (|55|l . then reduces the expressions (|58|l to 
those reported in the Appendix (apart from the overall 
factor — TO^a_F/(87r)). As a consequence, the expression 
for n' coincides within the two (Bogoliubov and Popov) 
approximations in strong coupling. 

The final form of the fermionic self-energy is eventually 
obtained by reconsidering the expressions H33|l . where the 
Popov propagator F^"^ {q) of Eq. (|5^ takes the place of 
Tii^q) and A satisfies the condition H55|l in the place of 
the gap equation 14911 . In this way, in the strong-coupling 



limit the composite bosons are described by the Popov 
approximation in the place of the Bogoliubov approxi- 
mation of Section II. 

The last step of our program would require us to mod- 
ify also the density equation (|34|l . consistently with the 
above approximations. For a homogeneous system, how- 
ever, no difference would show up in strong coupling from 
the result obtained in subsection IIC. 

All the above results within the Popov approxima- 
tion for composite bosons are consistent with the rela- 
tion as = 2ap. The same treatment which was made in 
subsection IIB to improve on this relation (in order to 
reduce the value of the ratio as/ap from 2 to 0.75) can 
be here applied. What one should do in this case is to 
first rewrite 'S^^{q)ii of Eq. 1)52(1 in terms of the (bare) 
boson-boson interaction H27|l as follows 



^Pop 



11 - ^ 2 / T^-^- i U2{q', q, q' , q) Fii(g') , 



(27r)3 13 



(59) 

and then replace U2 herein by the t-matrix for the 
composite bosons of Eq. I(31|l keeping the four- vector ar- 
guments of Eq. H59|l . The resulting bosonic self-energy 
then adds to the diagonal components of Eq. ((29|l when 
the corresponding replacement is made. The Popov prop- 
agator for the composite bosons is then obtained as in 
subsection IIB via the Dyson's equation, in a similar 
fashion to the procedure for point-like bosons described 
in the Appendix. Correspondingly, the gapless condition 
is given by Eq. I|30|l where now the diagonal self-energy 
includes also the Popov contribution with ts- This condi- 
tion reduces again to the gap equation H48|l in the strong- 
and weak-coupling limits (albeit with the modified value 
of as in the strong-coupling limit), in analogy to what 
was already shown for Eq. (|55() . 

It is worth emphasizing that, similarly to the anal- 
ogous treatment of subsection IIB within the Bogoli- 
ubov approximation, the price to pay for including the 
boson-boson scattering processes beyond the Born ap- 
proximation is to give up the self-consistency of the 
fermionic single-particle Green's functions entering the 
ladder propagators ((15|l . This self-consistency is actu- 
ally irrelevant in the strong-coupling limit, since it intro- 
duces subleading fermionic processes which do not map 
in that limit into corresponding processes for the com- 
posite bosons, in analogy to what was shown to occur in 
the normal phasei^. 



IV. CONCLUDING REMARKS 

In this paper, a diagrammatic theory for the BCS-BEC 
crossover has been set up for the constituent fermions, 
that treats the composite bosons which form in strong 
coupling in analogy to the Popov approximation for 
point-like bosons, with mass mg — 2m, chemical poten- 
tial = 2/x-|-eo, and coupling constant AiraB /tub- Care 
has also been taken to include diagrammatic terms which 
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modify the value of the boson-boson scattering length qb 
from the Born result as — '^o-Fi following the procedure 
introduced in Ref. 

A method of successive improvements of the fermionic 
diagrammatic theory has been presented, according to 
which the Dyson's equations for the constituent fermions 
and for the composite bosons have been treated at pro- 
gressive levels of sophistication, in an open-ended fashion. 
At any step, physical quantities (like the density) which 
can be calculated directly from the fermionic single- 
particle Green's functions then include more physical 
processes than the propagator for the composite bosons 
entering the fermionic self-energy. It is thus clear that the 
relation between conserving and gaplcss approximations 
for composite bosons, which was formally established in 
Ref. does not apply by construction to the present 
scheme of successive approximations. Nevertheless, it is 
significant that the bosonic self-energy of Fig. 4 results 
from taking partially into account the dressing of the 
fermionic single-particle lines within the self-consistent 
fermionic t-matrix approximation. It can also be shown 
that additional diagrams for the effective interaction in 
the particle-particle channel, which are required to im- 
plement the conserving nature of this fermionic t-matrix 
approximation in this channel^^i, can be associated with 
(at least some of) the boson-boson scattering processes 
beyond the Born approximation discussed in Ref. 1T3 . 

In the proposal presented in this paper, we have 
required the composite-boson propagator entering the 
fermionic self-energy in the broken-symmetry phase to 
be gapless. Such an approximation, however, cannot be 
conserving at the same time. The reason for this choice is 
that for the BCS-BEC crossover it is important to control 
the many-body approach in the two BCS and BEG limits, 
where the small parameter kp\aF\ is used to control the 
approximations. In this spirit, it appears appropriate to 
guarantee the composite-boson propagator entering the 
fermionic self-energy to be gapless at any coupling, while 
satisfying the conserving criterion only perturbatively in 
the small parameter. This approach has to be contrasted 
with the use of the self-consistent fermionic t-matrix ap- 
proximation in the broken-symmetry phase to deal with 
the BGS-BEG crossoveriS. In that approach, while sat- 
isfying exactly the conserving criterion at the single- 
particle level for any coupling, the composite-boson prop- 
agator entering the fermionic self-energy is gapped. This 
is because the additional diagrams for the effective in- 
teraction in the particle-particle channel, which would 
be needed to make this propagator gapless^^, are not 
included when constructing this self-consistent t-matrix 
self-energy. 

It is important to emphasize that the fermionic self- 
energy of Fig. 4, which was exploited in Section III 
to build up the Popov approximation for the compos- 
ite bosons from the fermionic two-particle Green's func- 
tions, is associated with pseudo-gap phenomena for the 
constituent fermions, both in the normal^ and broken- 
symmetry'*'^'* phases. In this way, the refinements of 



the Bogoliubov approximation for the composite bosons, 
which led to the Popov approximation, require one to in- 
clude pseudo-gap phenomena not only in the (fermionic) 
single-particle Green's function but also in the two- 
particle Green's function. In a related fashion, as pseudo- 
gap phenomena evolve with continuity across the critical 
temperature, the Popov theory for the composite bosons 
in the broken-symmetry phase similarly evolves into the 
theory of dilute interacting composite bosons in the nor- 
mal phasei4. 

The Popov approximation in the broken-symmetry 
phase is known to be valid^ for temperatures at which 
the condensate (uq) and noncondensate (n') densities are 
comparable to each other and to the total density n. This 
excludes the temperature ranges close to zero tempera- 
ture and to the critical temperature, where the Popov 
approximation may produce incorrect results. Further 
improvements beyond the Popov approximation are ac- 
cordingly required in these temperature ranges. 

The theoretical framework presented in this paper for 
the BCS-BEC crossover within the Popov approximation 
for the composite bosons awaits implementation by ex- 
plicit numerical calculations of physical quantities, simi- 
larly to what was recently done within the Bogoliubov ap- 
proximation for the composite bosons, both for a homo- 
geneous systemiiM and for trapped atomsii. At least in 
strong coupling, differences between the Popov and Bo- 
goliubov approximations should show up only in the pres- 
ence of a trapping potential, since for a homogeneous sys- 
tem these approximations differ only as far as the value of 
the bosonic chemical potential is concerned. Work along 
these lines is in progress. Numerical work should, in par- 
ticular, assess the deviations from the Popov approxima- 
tion for point-like bosons due to the composite nature of 
the bosons when approaching the intermediate-coupling 
region, an issue which can be subject to experimental 
verification. In this context, one may preliminary test to 
what an extent the occurrence of a downward shift in the 
critical temperature of point-like bosons treated within 
the Popov approximation**^ persists when the composite 
nature of the bosons begins to matter. 
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APPENDIX A: BOGOLIUBOV AND POPOV 
APPROXIMATIONS FOR POINT-LIKE BOSONS 

For the sake of completeness, in this Appendix we 
briefly recall the Bogoliubov and Popov approximations 
for point-like bosons in the homogeneous case. The re- 
sults here reported are used in the text to obtain simi- 
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lar descriptions for the composite bosons in the strong- 
coupHng limit of the fermionic attraction. 

Quite generally, the Dyson's equation for the propaga- 
tors Gb((7) of point-like bosons in the broken-symmetry 
phase reads (in matrix notation): 

Gsiq) = Gf{q) + Gf{q)-EB{q)GB{q) • (Al) 

In this expression, G^g\q) is the free-boson propagator 
with inverse 



goliubov approximation has been used and where 



Gfiq)-' 











(A2) 

where ms and fis are the bosonic mass and chemical 
potential, respectively, and ^^(g) is the 2x2 bosonic 
self-energy. 

Within the Bogoliubov approximation, S b (q) is taken 
of the formi 

^fHq) = 5^0 2) (^3) 

where rig is the condensate density. In this expression, 
the constant coupling g = AiraB/mB results from replac- 
ing the bare (repulsive) boson-boson potential by a t- 
matrix, which is in turn considered equal to the constant 
value for the scattering of two bosons in vacuum asso- 
ciated with the scattering length as^^- With the self- 
energy IjASp . the Dyson's equation (|A1|) yields: 



GfHq) 



1 



^l + EBiqr 



^ 2m B 



gn-Q 



(A4) 



Eb{ci) 




gn-o 



52 nl (A5) 



is the (gapless) Bogoliubov quasiparticle dispersion. 

The Popov approximation is obtained by considering 
the additional diagonal self-energy^ 



(g) = 25n'(^J 5) (A6) 



where the relation ^b — gn-o that holds within the Bo- 



where n' is the noncondensate density such that uq + n' 
is the total bosonic density. The bosonic propagators 
G^°'' (q) within the Popov approximation can be related 
to the propagators ljA4p within the Bogolibov approxi- 
mation via the modified Dyson's equation: 

Gl°^{q) = GfHq) + G|°s(,)sP°P(g)GP°P(<j) . (A7) 

The relation 11 b = gi^o + 2n') is required for the prop- 
agators to be gapless within the Popov approximation. 
Exploiting this relation, it can be shown that the Popov 
propagators G^°^(q) [Eq. ljA7p ] acquire the same form 
of the Bogohubov propagators G^°^ (g) [Eq. HA4|l ] , once 
the respective chemical potentials are eliminated in fa- 
vor of the condensate density uq. For the homogeneous 
case, the expression for n' thus coincides within the two 
(Bogoliubov and Popov) approximations'^^. 

The above coincidence between the Bogoliubov and 
Popov bosonic propagators stems from the simplifying 
assumption that the self-energies (|A3|I and l|A6p contain 
the constant coupling g, thus omitting the complicated 
(frequency and wave-vector) structure of the bosonic 
many-body t-matrix occurring in the formal expressions 
of the self-energies reported in Ref|5|. 
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